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^ (— I Abstract. We review the main features and results of thermal leptogenesis within 

f> | the type I seesaw mechanism, the minimal extension of the Standard Model explaining 

^ ^ neutrino masses and mixing. After presenting the simplest approach, the vanilla 

D , scenario, we discuss various important developments in recent years, such as the 

inclusion of lepton and heavy neutrino flavour effects, a description beyond a 
hierarchical heavy neutrino mass spectrum and an improved kinetic description 
within the density matrix and the closed-time-path formalisms. We also discuss how 
leptogenesis can ultimately represent an important phenomenological tool to test the 
seesaw mechanism and the underlying model of new physics. 
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1. Introduction: the double side of leptogenesis 

A successful model of baryogenesis cannot be realised within the Standard Model 
(SM) and, therefore, the observed matter-antimatter asymmetry of the Universe can 
be regarded as an evidence for new physics beyond the SM. 

The discovery of neutrino masses and mixing in neutrino oscillation experiments in 
1998 pQ has for the first time shown directly, in particle physics experiments, that the 
SM is indeed incomplete, since it strictly predicts, in the limit of infinite cutoff, that 
neutrinos are massless and, therefore, cannot oscillate. 

This discovery has greatly raised the interest for leptogenesis [21 E], a model of 
baryogenesis that is a cosmological consequence of the most popular way to extend the 
SM in order to explain why neutrinos are massive but at the same time much lighter 
than all the other fermions: the seesaw mechanism [4]. As a matter of fact, leptogenesis 
realises a highly non-trivial link between two completely independent experimental 
observations: the absence of primordial antimatter in the observable Universe and the 
observation that neutrinos mix and (therefore) have masses. In fact, leptogenesis has a 
naturally built-in double-sided nature. On one side, it describes a very early stage in 
the history of the Universe characterised by temperatures T\ ep > 100 GeV, much higher 
than those probed by Big Bang Nucleosynthesis, T B bn ~ (0-1-1) MeV; on the other side, 
it complements low-energy neutrino experiments providing a completely independent 
phenomenological tool to test models of new physics embedding the seesaw mechanism. 

In this article we review the main features and results of leptogenesis. Let 
us give a brief outline. In Section 2 we present the status of low-energy neutrino 
experiments measuring neutrino masses and mixing parameters, and we introduce the 
seesaw mechanism, which provides an elegant framework to explain them. In Section 
3 we discuss the vanilla leptogenesis scenario. In Section 4 we show the importance 
of accounting for flavour effects for a correct calculation of the final asymmetry. In 
Section 5 we discuss the density matrix formalism which properly takes into account 
decoherence effects, which are crucial to describe the transition from a one-flavoured 
regime to a fully flavoured regime. In Section 6 we relax the assumption of hierarchical 
right-handed neutrino mass spectrum, and discuss how the asymmetry can be calculated 
in the degenerate limit. In Section 7 we discuss different ways to improve the kinetic 
description beyond the density matrix formalism. In Section 8 we discuss other effects 
(thermal corrections, spectator processes, scatterings) that have been considered and 
that can give in some cases important corrections. In Section 9 we show how leptogenesis 
represents an important guidance to test models of new physics. Finally, in Section 10 
we conclude outlining the prospects to test leptogenesis in future years. 

2. Neutrino masses and mixing 

Neutrino oscillation experiments have established two fundamental properties of 
neutrinos. The first one is that the neutrinos mix. This means that the neutrino 



The minimal scenario of leptogenesis 



3 



weak eigenfields v a [a = e, //, r) do not coincide with the neutrino mass eigenfields z/j 
(i = 1,2,3) but are obtained applying to them a unitary transformation described by 
the (3 x 3) leptonic mixing matrix U, 

i 

The leptonic mixing matrix is usually parameterised in terms of 6 physical parameters, 
three mixing angles, 9\ 2) #13 and #23, and three phases, two Majorana phases, p and a 
and one Dirac phase S, 

^ C12 C13 S 12 C13 Si3e lS \ 

U = 



-S12C23 - C 12 S 2 3S 13 e l5 C 12 C 23 - S 12 S 23 S 13 e lS S23C13 

y S12 S23 — C12 C23 S13 e 1<5 — c 12 s 2 3 — s 12 C23 S13 e lS C23C13 J 



diag(e^,l,e^) ,(2) 



where = sin 6^- and Cy = cos%. A global analysis [5] of all existing neutrino data, 
prior to the results of a non- vanishing # 13 from short baseline reactors, gives Q\ 2 = 34°±1° 
for the solar mixing angle, #23 = 40.4°^g O for the atmospheric mixing angle, and 
#13 = 9.0° ± 1.3° for the reactor mixing angle, where the latter is mainly dominated 
by the evidence for a non- vanishing # 13 found by the T2K experiment [6] that confirmed 
previous hints {?]. Recently, the Daya Bay and Reno short-baseline reactor neutrino 
experiments respectively found 6 13 = 8.8° ± 0.8° ± 0.3° [BJ, and 13 = 9.8° ± 0.6° ± 0.85° 
[9], confirming, at more than 5a, the previous results. 

The second important property established by neutrino oscillation experiments is 
that neutrinos are massive. More specifically, defining the three neutrino masses in a 
way that mi < m 2 < 777,3, neutrino oscillation experiments measure two mass-squared 
differences that we can indicate with Am^ m and Am s 2 oI since historically the first one 
has been first measured in atmospheric neutrino experiments and the second one in solar 
neutrino experiments. Two options are currently allowed by previous experiments. A 
first option is 'normal ordering' (NO) and in this case 

m 2 3 - m 2 2 = Am 2 atm and m\-m\ = Am s 2 ol , (3) 

while a second option is represented by 'inverted ordering' (IO) and in this case 

m\ - m\ = Am 2 sol and m\-m\ = Am^ tm . (4) 

It is convenient to introduce the atmospheric neutrino mass scale m atm = 
a/ Am atm + Am 2 ol = (0.049 ± 0.001) eV and the solar neutrino mass scale m so \ = 
V'AmL = (0.0087 ± 0.0001) eV [5]. 

The measurements of m atm and m so \ are not sufficient to fix all three neutrino 
masses. If we express them in terms of the lightest neutrino mass m\ we can see from 
Fig. 1 that while 777.2 > "^soi and m 3 > m atm , the lightest neutrino mass can be arbitrarily 
small implying that the lightest neutrino could be even massless. 

The lower limits for m 2 and m 3 are saturated when m\ m so \. In this case one has 
hierarchical neutrino models, either normal and in this case m 2 ~ m so \ and m 3 ~ ?n atm , 
or inverted, and in this case m 2 ~ a/ m 2 tm — 777g 0l ~ m 3 ~ m atm . On the other hand, for 
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Figure 1. Neutrino masses rm versus the lightest neutrino mass mi. The three 
upper bounds (at 95% C.L.) discussed in the body text from absolute neutrino 
mass scales phenomenologies are also indicated. 

mi 3> m atm one obtains the limit of quasi-degenerate neutrinos when all three masses 
can be arbitrarily close to each other. 

However, the lightest neutrino mass is upper bounded by absolute neutrino mass 
scale experiments. Tritium beta decay experiments [10J place an upper bound on the 
effective electron neutrino mass m Ue < 2eV (95% C.L.) that translates into the same 
upper bound on ni\. This is derived from model- independent kinematic considerations 
that apply independently of whether neutrinos have a Dirac or Majorana nature. 

Neutrinoless double beta decay (0v{3{3) experiments place a more stringent upper 
bound on the effective Qvf3f3 Majorana neutrino mass, m ee < (0.34-0.78) eV (95% C.L.) 
as obtained by the CUORICINO experiment [TT] 0. It translates into the following 
upper bound on mi [T2] : 

mi < m ee /(cos 26 12 cos 2 9 13 - sin 2 13 ) < 3.45 m ee < (1.2-2.7) eV (95% C.L.) . (5) 

Here, the wide range is due to theoretical uncertainties in the calculation of the involved 
nuclear matrix elements. However, this upper bound applies only if neutrinos are of 
Majorana nature, which is the relevant case for us, since the seesaw mechanism predicts 
Majorana neutrinos. 

1 At 90% C.L. the CUORICINO result is m ee < (0.27-0.57) eV. For comparison, the bound from the 
Heidelberg-Moscow experiment is m ee < (0.21-0.53) eV (90% C.L.) [13] while the EXO-200 experiment 
has recently found the upper bound m ee < (0.14-0.38) eV (90% C.L.) [M]. 
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From cosmological observations, within the ACDM model, one obtains a very 
stringent upper bound on the sum of the neutrino masses [15] 

^mi<0.58eV (95%C.L.). (6) 

i 

This translates into the most stringent upper bound that we currently have on the 
lightest neutrino mass m x < 0.19 eV (95% C.L.), an upper bound that almost excludes 
quasi-degenerate neutrino models §. 

A minimal extension of the SM, able to explain not only why neutrinos are massive 
but also why they are much lighter than all the other massive fermions, is represented 
by the seesaw mechanism [4j. In the minimal type I version, one adds right-handed 
neutrinos Nm to the SM lagrangian with Yukawa couplings h and a Majorana mass 
term that violates lepton number 

£ = C SM + i N^i^NiR - J^l Ki N iR $--M i Nf R N iR + h.c. (7) 

For definiteness we will consider the case of three RH neutrinos (i = 1,2,3). This is 
the most attractive case corresponding to have one RH neutrino for each generation, 
as for example nicely predicted by 5*0(10) grand-unified models. Notice, however, that 
all current data from low-energy neutrino experiments are also consistent with a more 
minimal model with only two RH neutrinos that will be discussed in Section 9. 

After spontaneous symmetry breaking, a Dirac mass term mjj = v h is generated 
by the Higgs vev v. In the seesaw limit, M 3> mo, the spectrum of neutrino masses 
splits into a light set given by the eigenvalues mi < vri2 < of the neutrino mass 
matrix 

m v = -m D ^- m T D , (8) 

and into a heavy set M\ < M 2 < M 3 coinciding to a good approximation with the 
eigenvalues of the Majorana mass matrix corresponding to eigenstates Ni ~ N iR + N? R . 
The symmetric neutrino mass matrix m v is diagonalized by a unitary matrix U, 

D m = di&g(m 1 ,m 2 ,m 3 ) = -Ulm u U*. (9) 

In a basis where the charged lepton mass matrix is diagonal, U v coincides with the 
leptonic mixing matrix U. In this way the lightness of ordinary neutrinos is explained 
just as an algebraic by-product. If the largest eigenvalue in the Dirac neutrino mass 
matrix is assumed to be of the order of the electroweak scale, as for the other massive 
fermions, then the atmospheric neutrino mass scale m a t m can be naturally reproduced 
for M3 ~ 10 14-15 GeV, very close to the grand-unified scale. This is the minimal version 
of the seesaw mechanism, often indicated as type I seesaw mechanism. 

The seesaw formula (jS]) can be recast as an orthogonality condition for a matrix Vt 
that in a basis where simultaneously the charged lepton mass matrix and the Majorana 

2 For a more general discussion of neutrino mass bounds, in particular on theoretical assumptions and 
uncertainties, see for example [T2 l HB j ITT]. 
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mass matrix are diagonal, provides a useful parametrisation of the neutrino Dirac mass 
matrix [18J, 



where D M = diag(M 1; M 2 , M 3 ). The orthogonal matrix contains 6 independent 
parameters and together with the neutrino masses, determines the properties of the 
heavy RH neutrinos, the three lifetimes and the three total CP asymmetries. In this 
way, one can easily see that the total number of additional parameters introduced by the 
see-saw lagrangian is 18 in the case of 3 RH neutrinos. The orthogonal parametrisation is 
quite useful since, given a model that specifies and the three RH neutrino masses Mj, 
one can easily impose both the experimental information on the 9 low-energy neutrino 
parameters (6 contained in U and the 3 mj) and, as we will see, also the requirement of 
successful leptogenesis. 

3. The simplest scenario: vanilla leptogenesis 

Leptogenesis belongs to a class of models of Baryogenesis where the asymmetry is 
generated from the the out-of-equilibrium decays of very heavy particles. Interestingly, 
this is the same class to which the first model proposed by Sakharov belongs. This 
class of models became very popular with the advent of Grand Unified Theories 
(GUT) that provided a well-defined and motivated framework. In GUT baryogenesis 
models the very heavy decaying particles generating the asymmetry are the same new 
gauge bosons predicted by the GUT's. However, the final asymmetry depends on 
too many untestable parameters, so that imposing successful baryogenesis does not 
lead to compelling experimental predictions. This lack of predictability is made even 
stronger considering that the decaying particles are too heavy to be produced thermally 
and one has, therefore, to invoke a non-thermal production mechanism of the gauge 
bosons. This is because while the mass of the gauge bosons is about the GUT scale, 
Mx ~ 10 15-16 GeV, the reheating temperature at the end of inflation T re h cannot be 
higher than ~ 10 15 GeV from CMB observations, where the reheating temperature is 
the value of the temperature at the beginning of the radiation-dominated regime after 
inflation [T9] . 

The minimal (and original) version of leptogenesis, that we discuss in this review, 
is based on the type I seesaw mechanism [2], where the asymmetry is produced by the 
three heavy RH neutrinos predicted by the see-saw mechanism and whose masses can 
have values orders-of-magnitude below the upper bound on the T rc h < 10 15 GeV from 
inflation and CMB observations. We will call 'minimal leptogenesis scenarios' those 
scenarios based on a type I seesaw mechanism and on a thermal production of the 
RH neutrinos, implying that T re h cannot be too much below (at least) the lightest RH 



neutrino mass M\ q At these high temperatures the RH neutrinos can be produced by 

3 For a description of leptogenesis from RH neutrino oscillations, and for leptogenesis scenarios based 
on models beyond the type I see-saw mechanism, see respectively [20] and [21] in this Focus Issue. 



m D = UD 




(10) 
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the Yukawa interactions of leptons and Higgs bosons in the thermal bath. A solution 
of the Boltzmann equations shows that, more exactly, T rc h can be even up to 10 times 
lower than Mi, the exact value depending on the strength of the coupling [22] . 

After their production, the RH neutrinos decay either into lepton and Higgs 
doublets Ni — > £i + $ (when leptogenesis takes place, the electroweak symmetry is 
not broken) with a decay rate T, or into anti-lepton and conjugate Higgs doublets 
Ni — > ii + $^ with a decay rate F. Both inverse processes and decays violate lepton 
number (AL = 1) and in general CP as well. It is also important to notice that they 
also violate the B — L asymmetry, (A(5 — L) = 1). At temperatures T ^> 100 GeV, 
non-perturbative Standard Model processes called sphalerons are in equilibrium [23J. 
They violate both lepton and baryon number while they still conserve B — L. Hence 
the lepton asymmetry produced in the elementary processes is reprocessed in a way 
that ~ 1/3 of the B — L asymmetry is in the form of baryon number while ~ —2/3 
of the B — L asymmetry is in the form of a lepton number [24J. Therefore, two of the 
Sakharov conditions are satisfied. The third Sakharov condition, departure from thermal 
equilibrium, is also satisfied since a fraction of the decays occur out of equilibrium, 
implying that part of the asymmetry survives the washout from inverse processes. 

In the most general approach, even within minimal leptogenesis, the asymmetry 
depends on all the 18 seesaw parameters and the calculation itself presents different 
technical difficulties. However, there is a simplified scenario p[2^[2^[271|28l[29|[22|[3n] 
that grasps most of the main features of leptogenesis and is able to highlight important 
connections with the low-energy neutrino parameters in an approximated way. We will 
refer to this scenario as 'vanilla leptogenesis'. Let us discuss the main assumptions and 
approximations. 

The first assumption is the one-flavour approximation, where the flavour 
composition of the leptons produced by (or producing) the RH neutrinos has no influence 
on the final value of the asymmetry and can be therefore neglected. It also assumes that 
leptons produced by different RH neutrinos can be treated, for all practical purposes, 
as having the same flavour (i.e. t\ = £2 = £3 = £)■ This is equivalent to saying that 
what counts is just the asymmetry between the total number of leptons and the total 
number of anti-leptons irrespectively of whether these leptons are electron, muon or 
tauon doublets, and irrespectively of the RH neutrino that generated it. Under this 
assumption the final B — L asymmetry is the sum of two contributions 

NL L = Kr L ^ + N^j L . (ii) 

The first term on the right-hand side, A^ r _!^ ex,f , is the residual value of a possible 
pre-existing asymmetry generated by some external mechanism prior to the onset of 
leptogenesis. It has to be regarded as a possible external contribution whose initial value 
would be set up by some independent source of leptogenesis. Therefore, if leptogenesis 
is responsible for the observed asymmetry, this term has to be negligible. 

The second term is the genuine leptogenesis contribution to the final asymmetry 



The minimal scenario of leptogenesis 



8 



and is given by the sum of the three contributions from each RH neutrino species, 

N^J L = Y,^4(Ki,K 2} K 3 ). (12) 

i 

Each contribution is the product of the total CP asymmetry Si times the final value of 
the efficiency factor, «f(.Ki, K 2 , K 3 ), which depends on the decay parameters K iy defined 
as Ki = Fx)i(T = 0)/H(T = Mj), where H is the Hubble expansion rate, and the total 
decay rates are Td^ = Tj + f j. The total CP asymmetries are defined as 

(13) 

and a perturbative calculation from the interference of tree level with one loop self- 
energy and vertex diagrams gives [HH E2J [33] 



_3_ ^ Im [(/it ] 

3T* 



h)u ^xj/xi ' 



having introduced 

ax) 



' X 



'1 + x) In 



1 + x 



X 



X 



1 — X 



(14) 



(15) 



and defined Xj = M^jM\. The efficiency factors are computed solving a simple set of 
Boltzmann equations integrated over the momenta (rate equations) 
dN Ni 



dz 
dN B _ L 
dz 



A (N Ni - N%) 



1,2,3 



i=i 



AW(*) + ]T w! D ( 



(16) 
(17) 



where 2; = Mi/T and where we indicated with Nx any particle number or asymmetry 
X calculated in a portion of co-moving volume containing one heavy neutrino in ultra- 
relativistic thermal equilibrium, so that N^.(T 3> Mj) = 1. With this convention, the 
predicted baryon-to-photon ratio tjb is related to the final value of the final B — L 
asymmetry by the relation 



Vb = asph^^ ^ 0.96 x l(T 2 iV f 

7 



B-L 5 



'II 



where N* ec ~ 37, and a sp h = 28/79. Defining Z{ = z yfxl, the decay factors are defined 
as 



D, 



' L> "' K-iXiZ { — 



^ , , (19) 

Hz \ 74 / 

where the thermally averaged dilation factors (I/7) are given by the ratio K,i(z) / K, 2 (z) 
of the modified Bessel functions. After proper subtraction of the resonant contribution 
from AL = 2 processes in order to avoid double counting [341 129] , the inverse decay 
washout terms are simply given by 

1 



W} D (z 



KiJxltCAzAz; 



(20) 
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The washout term AW(z) is the non- resonant AL = 2 processes contribution. It gives 
a non-negligible effect only at z ^> 1 and in this case it can be approximated as [27] 



z 2 Vl0 10 GeVy VeV 2 ^ 

where u> ~ 0.186 and m 2 = m\+m\+m\. The efficiency factors can be simply calculated 
analytically as 

K*(K t ) = - r dz'^e-f~ W'HE. (22) 

The second assumption is that the RH neutrino mass spectrum is hierarchical with 
M 2 >3Ml 

Under these two assumptions (i.e. unflavoured description plus hierarchical RH 
neutrino mass spectrum) and barring a particular case, that we will discuss later on, 
the final asymmetry is typically dominantly produced by the lightest RH neutrino out- 
of-equilibrium decays, in a way that the sum in Eq. ( fTTT) can be approximated by the 
first term (i = 1), explicitly 

JVL. L =* 6! (23) 

so that an 'iVj-dominated scenario' is realised. This happens either because \£2,3\ "C 
or because the asymmetry initially produced by the N 2t 3 decays is afterwards washed out 
by the lightest RH neutrino inverse processes such that k 23 (K 2 ,3) "C k\(K\). Indeed, 

(2 3*1 

if we indicate with N B !_[(T > Mi) the contribution to the Nb-l asymmetry from the 
two heavier RH neutrinos prior to the lightest RH neutrino washout, the final values 
are given simply by 

NtY = Nt 3) L (T>M 1 )e-% K i. (24) 

The same exponential washout factor also suppresses the residual value of a possible pre- 
existing asymmetry. Therefore, it is sufficient to impose the strong washout condition 
Ki 3> 1 for the pre-existing asymmetry and for the contribution from heavier RH 
neutrinos to be negligible. The decay parameters can be expressed in terms of the 
see-saw parameters as 

(25) 

where ~ 1.1 x 10 -3 eV is the equilibrium neutrino mass EH]- It is therefore 
quite straightforward that, barring cancellations in the seesaw formula, one typically 
has Ki ~ i^soi-i^atm ~ 10-50 > 1, where if so i( a tm) = ^soi(atm)/^- The same condition 
also guarantees that the final asymmetry is independent of the initial 7Vi-abundance. It 
is then quite suggestive that the measured values of m sol and m atm have just the right 
values to produce a washout that is strong enough to guarantee independence of the 
initial conditions but still not that strong to prevent successful leptogenesis [28]. This 
leptogenesis conspiracy between experimental results and theoretical predictions is the 
main reason that has determined the success of leptogenesis during the last years. 
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Figure 2. Left: Neutrino mass bounds in the vanilla scenario. Right: Relaxation 
of the lower bound on Mi thanks to additional unbounded flavoured CP violating 
terms [30] . 



There is a particular case where K\ ^> 1 does not hold and in this case the final 
asymmetry can be dominated by the contribution coming from the next-to-lightest RH 
neutrinos [35]. For this case one still has K% ^> 1, so that the independence of the initial 
conditions holds nonetheless. For the time being, as an additional third assumption, we 
will not consider this particular case. 

If one excludes cancellations among the different terms contributing to the neutrino 
masses in the seesaw formula, one obtains the Davidson-Ibarra upper bound on the 
lightest RH neutrino CP asymmetry [36] 

< max in-6 M l m atm , . 

£l - Sl " iU 10"> GeV m, + m 3 " 1 ] 

Imposing rf^ ~ 0.01 e™ 3 * k\ > ^§ MB , one obtains the allowed region in the plane 
(mi, Mi) shown in the left panel of Fig. 2. One can notice the existence of an upper 
bound on the light neutrino masses m\ < 0.12 eV [37:, 28J, incompatible with quasi- 
degenerate neutrino mass models, and a lower bound on M\ >3x 10 9 GeV [361 EH I2H] 
implying a lower bound on the reheat temperature T re h > 10 9 GeV [22J. 

An important feature of vanilla leptogenesis is that the final asymmetry does not 
directly depend on the parameters of leptonic mixing matrix U. That implies that one 
cannot establish a direct model-independent connection. In particular the discovery of 
CP violation in neutrino mixing would not be a smoking gun for leptogenesis and, vice 
versa, a non-discovery would not rule out leptogenesis. However, within more restricted 
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scenarios, for example imposing some conditions on the neutrino Dirac mass matrix, 
links can emerge. In Section 9 we will discuss in detail the interesting case of 5*0(10)- 
inspired models. 

In the past years different directions beyond the vanilla leptogenesis scenario were 
explored, usually aiming at evading the above-mentioned bounds. For example, it was 
noticed that the Davidson- Ibarra bound eq. fl26l) is strictly speaking evaded by an extra- 
term contribution As\ to the total CP asymmetry [381 EQ] ■ However, this extra-term 
vanishes if M3 ~ M2 and is suppressed as Ae\ oc (M1/M2) 2 . It can significantly relax 
the bounds on neutrino masses only when |f2jj| 2 > (M2/M1) 2 , implying cancellations 
with a certain degree of fine tuning in the see-saw formula for the neutrino masses. 
For example, in usual models with |f^j| 2 < 10, this extra-contribution can be safely 
neglected for M 2 > 10 Mi. 

The development that proved to have the most important impact on the final 
asymmetry, compared to a calculation within the vanilla scenario, is certainly the 
inclusion of flavour effects. For this reason we discuss them in detail in the next Section. 

4. The importance of flavour 

The inclusion of flavour effects provides the most significant modification in the 
calculation of the final asymmetry compared to the vanilla scenario. Two kinds of 
flavour effects are neglected in the vanilla scenario: heavy neutrino flavour effects, how 
heavier RH neutrinos can contribute to the final asymmetry, and lepton flavour effects, 
how the flavour composition of the leptons quantum states produced in the RH neutrino 
decays affects the calculation of the final asymmetry. We first discuss the two effects 
separately, and then we show how their interplay can have very interesting consequences. 

4.1. Heavy neutrino flavour effects 

In the vanilla scenario the contribution to the final asymmetry from the heavier 
RH neutrinos is negligible because either the CP asymmetries are suppressed in the 
hierarchical limit compared to £™ ax (cf. Eq. (126)) ) and/or because, even assuming that a 
sizeable asymmetry (compared to the observed one) is produced at T ~ ^2,3, it is later 
on washed out by the lightest RH neutrino inverse processes. 

However, as we anticipated, there is a particular case when, even neglecting the 
lepton flavour composition and assuming a hierarchical heavy neutrino mass spectrum, 
the contribution to the final asymmetry from next-to-lightest RH neutrino decays 
can be dominant and explain the observed asymmetry [35]. This case corresponds 
to a particular choice of the orthogonal matrix such that Aq is so weakly coupled 
(corresponding to K\ <C 1) that its washout can be neglected. For the same choice 
of the parameters, the N2 total CP asymmetry 62 is unsuppressed if M3 < 10 15 GeV. In 
this case a iV 2 -dominated scenario is realised. Notice that in this case the existence of 
a third (heavier) RH neutrino species is crucial in order to have a sizeable 62- 
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The contribution from the two heavier RH neutrino species is also important in 
the quasi-degenerate limit when 5i = (Mj — Mi) /Mi < 1, i — 2,3. In this case the 
CP asymmetries £2,3 are not suppressed, and the washout from the lighter RH neutrino 
species is moderate, with no exponential prefactor [39| HO]. 



4-2. Lepton flavour effects 

The importance and generality of flavour effects in leptogenesis was fully highlighted 
in |JTJ H2]. Their role was first discussed in [26J and included in specific scenarios in 

pi mug. 

For the time being, let us continue to assume that the final asymmetry is dominantly 
produced from the decays of the lightest RH neutrinos Ni, neglecting the contribution 
from the decays of the heavier RH neutrinos N 2 and N 3 . If Mi ^> 10 12 GeV, the flavour 
composition of the quantum states of the leptons produced from Ni decays has no 
influence on the final asymmetry and a one-flavour regime holds [261 EJ E2j- This is 
because the lepton quantum states evolve coherently between the production from a 
iVi-decay and a subsequent inverse decay with a Higgs boson. In this way the lepton 
flavour composition does not play any role. 

However, if 10 12 GeV > Mi > 10 9 GeV, during the relevant period of generation of 
the asymmetry, the produced lepton quantum states will, on average, have an interaction 
with RH tauons before undergoing the subsequent inverse decay. In this way the tauon 
component of the lepton quantum states is measured by the thermal bath and the 
coherent evolution breaks down [26 | I¥T | H2]. Therefore, at the subsequent inverse decays, 
the lepton quantum states are an incoherent mixture of a tauon component and of a 
(still coherent) superposition of an electron and a muon component that we can indicate 
with t- 1 . 

The fraction of asymmetry stored in each flavour component is not proportional in 
general to the branching ratio of that component. This implies that the two flavour 
asymmetries, the tauon and the t 1 - components, evolve differently and have to be 
calculated separately. In this way the resulting final asymmetry can considerably differ 
from the result in the one-flavour regime. This can be indeed approximated by the 
expression 

Ni_ L ~ 2 6l k{Ki) + ^ [k(K w ) - k(K 1t )] , (27) 

where Ki a = p\ a K i: the p? a 's (a = r, r- 1 ) are the tree level probabilities that the leptons 
ii and the anti-leptons £[ produced in the decays of the Aq's are in the flavour a, while 
Api T is the difference between the probability to find ii in the flavour r and that one 
to find £'1 in the flavour f. If we compare this expression with the one-flavour regime 
result eq. ff23l one can see that if Ap lT = (or if K lT = K iT ±), then the final symmetry 
is enhanced just by a factor 2. However, in general leptons and anti-leptons have a 
different flavour composition [261 H2] and in this case Api T 7^ and the final asymmetry 
can be much higher than in the one-flavour regime. The most extreme case is when the 
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total B — L number is conserved (i.e. E\ = 0) and still the second term in eq. (I27p can 
be non- vanishing [32] and even explain the observed asymmetry [50J. 

If Mi < 10 9 GeV, even the coherence of the r 1 - component is broken by the muon 
interactions between decays and inverse decays and a full three flavour regime applies. 
In the intermediate regimes a density matrix formalism is necessary to properly describe 
decoherence [^T | 1361 I4T1 138] . 

We can briefly say that lepton flavour effects induce three major consequences that 
are all encoded in the expression ( 1271) valid in the two fully flavoured regime, i) The 
washout can be considerably lower than in the unflavoured regime [3D 132] since one 
can have that the asymmetry is dominantly produced in a flavour a with K\ a <C K\ 
ii) The leptonic mixing matrix enters directly in the calculation of the final asymmetry, 
more specifically and in particular the low-energy phases contribute as a second source 
of CP violation in the flavoured CP asymmetries [4*21 l4*9l [50j EI]- As an interesting 
phenomenological consequence, the same source of CP violation that could take place 
in neutrino oscillations could be sufficient to explain the observed asymmetry [501 I5"T] . 
though under quite stringent conditions on the RH neutrino mass spectrum and with 
an exact determination requiring a density matrix calculation [52]. Notice that this is 
a particular case realising the above-mentioned scenario with si = [42J. This problem 
becomes particularly interesting in the light of the recent discovery of a non- vanishing 
#13 angle [7J El M, E] , a necessary condition to have CP violation in neutrino oscillations. 
In particular a calculation of the lower bound on # 13 that necessarily requires the use of 
density matrix equations, iii) The flavoured CP asymmetries, given by 

- TE^m: £ { Im K<M^)«] ^Jg +35^1) to [KMh^] } . (28) 

contain a second term that conserves the total lepton number (it cancels exactly when 
summing over flavour in the total CP asymmetry and it, therefore, contribute to Api T in 
eq. ( l2~Tj) but not to £i), and therefore the upper bound in Eq. (}2l)|) does not strictly apply 
to the flavoured CP asymmetries. As a consequence, allowing for a mild cancelation 
in the neutrino mass matrix, corresponding to ~ 1, and also for a mild RH 

neutrino mass hierarchy (M2/M1 ~ 10), the lower bound on T re h can be relaxed by 
about one order of magnitude, down to 10 8 GeV [30], as shown in the right panel of 
Fig. 2. However, for many models such as sequential dominated models [53], these 
cancellations do not occur, and lepton flavour effects cannot relax the lower bound 
on T re h |50j. One known exception is given by the inverse seesaw model [53] which 
naturally explains large cancellations in the neutrino mass matrix and leads to large 
CP asymmetries thanks to an underlying lepton number symmetry. This leads to the 
relaxation of the lower bound on T reh by up to three orders of magnitude [551 ES] (see 
more detailed discussion in the next Section). 

Whether or not the upper bound m; < 0.1 eV on neutrino masses found in the 
vanilla scenario still holds in a flavoured iVi-dominated scenario, is relaxed, or even 
completely evaporates, is a controversial topic. It still surely holds in the one-flavour 
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Figure 3. The ten RH neutrino mass patterns corresponding to 
leptogenesis scenarios with different sets of classical Boltzmann equations 
for the calculation of the final asymmetry |58j . 



regime for Mi > 10 12 GeV but it certainly does not hold in the two-fully flavoured 
regime [UJ E7J [30]. However, it was found that the two-fully flavoured regime is not 
respected at large values mi > 0.1 eV [4"B] . In [U] it was found that it holds up to 
mi ~ 2eV, implying in any case an upper bound much above current experimental 
bounds and, therefore, uninteresting. In [30], including the information from low energy 
neutrino data and accounting for the Higgs asymmetry, it was found again that it holds 
only up to mi > 0.1 eV and in this case a density matrix approach would be required 
for a conclusion. 

4-3. The interplay between lepton and heavy neutrino flavour effects 

As we have seen, when lepton flavour effects are neglected, the possibility that the next- 
to-lightest RH neutrino decays contribute to the final asymmetry relies on a special 
case realising the iV 2 -dominated scenario [35J. On the other hand, when lepton flavour 
effects are taken into account, the contribution from heavier RH neutrinos cannot be 
neglected in a much more general situation. Even the contribution from the heaviest 
RH neutrinos can be sizeable (i.e. explain the observed asymmetry) and has to be taken 
into account in general. 

As a result, the calculation of the final asymmetry becomes much more involved. 
Assuming hierarchical mass patterns and that the RH neutrino processes occur in one 
of the three different fully flavoured regimes, one has to consider ten different mass 
patterns, shown in Fig. 2, that require specific multi-stage sets of classical Boltzmann 
equations for the calculation of the final asymmetry. 

4-3.1. The (flavoured) N 2 -dominated scenario Among these 10 RH neutrino mass 
patterns, for those three where the N\ washout occurs in the three-fully-flavoured regime, 
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for Mi <C 10 9 GeV, the final asymmetry has necessarily to be produced by the N 2 either 
in the one or in the two-fully-flavoured regime (for M 2 3> 10 9 GeV). They have some 
particularly attractive features and realise a 'flavoured ^-dominated scenario' jl5] (it 
corresponds to the fifth, sixth and seventh mass pattern in Fig. 3). 

While in the unflavoured approximation the lightest RH neutrino washout yields a 
global exponential washout factor, when lepton flavour effects are taken into account, the 
asymmetry produced by the heavier RH neutrinos, at the N\ washout, gets distributed 
into an incoherent mixture of charged lepton flavour eigenstates [35]- It turns out that 
the N\ washout in one of the three flavours is negligible, corresponding to have at least 
one Ki a < 1, in quite a wide region of the parameter space [30]. In this way, accounting 
for flavour effects, the region of applicability of the iV 2 -dominated scenario enlarges 
considerably, since it is not necessary that Ni fully decouples but it is sufficient that it 
decouples just in a specific lepton flavour |45j. The unflavoured ^-dominated scenario 
is recovered in the limit where all three Ki a < 1 and either M 2 > 10 12 GeV or for 
10 12 GeV > M 2 > 10 9 GeV and K 2 < 1. 

Recently, it has also been realised that, accounting for the Higgs and for the quark 
asymmetries, the dynamics of the flavour asymmetries couple and the lightest RH 
neutrino washout in a particular flavour can be circumvented even when Ni is strongly 
coupled in that flavour [59]. Another interesting effect arising in the iV 2 -dominated 
scenario is phantom leptogenesis. This is a pure quantum-mechanical effect that for 
example allows parts of the electron and of the muon asymmetries, the phantom terms, 
to undergo a weaker washout at the production than the total asymmetry. It has 
been recently shown that phantom terms associated to a RH neutrino species Ni with 
Mj ^> 10 9 GeV are present not just in the ^-dominated scenario [60] ■ However, it 
should be noticed that phantom terms produced by the lightest RH neutrinos cancel with 
each other and thus do not contribute to the final asymmetry, though they can induce 
flavoured asymmetries much larger than the total asymmetry, something potentially 
relevant in active-sterile neutrino oscillations [61J. 

4-3.2. Heavy neutrino flavour projection Even assuming a strong RH neutrino mass 
hierarchy, a coupled N% in all lepton flavours {K\ a ^> 1 for any a) and Mi > 10 12 GeV 
(it corresponds to the first panel in Fig. 2), the asymmetry produced by the heavier RH 
neutrino decays at T ~ Mi, in particular by the N 2 s decays, can be large enough to 
explain the observed asymmetry by avoiding most of the washout from the lightest RH 
neutrino processes. This is because, in general, there is an 'orthogonal' component that 
escapes the Ni washout [261 162] while the remaining 'parallel' component undergoes the 
usual exponential washout. For a mild mass hierarchy, 5 3 < 10, even the asymmetry 
produced by the N^s decays can be large enough to explain the observed asymmetry 
and escape the Ni and N 2 washout. Heavy neutrino flavour projection is also occurring 
when 10 12 GeV > Mi > 10 9 GeV, in this case in the e-fi plane. 

When the effect of heavy neutrino flavour projection is taken into account jointly 
with an additional contribution to the flavoured CP asymmetries e 2a that is not 
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suppressed when the heaviest RH neutrino mass M 3 > 10 15 GeV, this can lead to the 
possibility of a dominant contribution from the next-to-lightest RH neutrinos even in 
an effective two RH neutrino model that can be regarded clS cLS cL limit case in a (more 
appealing) three RH neutrino case with M 3 > 10 15 GeV [59J. 

4-3.3. The problem of the initial conditions in flavoured leptogenesis As we have 
seen, in the vanilla scenario the unflavoured assumption reduces the problem of the 
dependence on the initial conditions to simply imposing the strong washout condition 
Ki 1. In other words, there is a full equivalence between strong washout and 
independence of the initial conditions. 

When (lepton and heavy neutrino) flavour effects are considered, the situation is 
much more involved and for example imposing strong washout conditions on all flavoured 
asymmetries (K ia ) is not enough to guarantee independence of the initial conditions. 
Perhaps the most striking consequence is that in a traditional iVi-dominated scenario 
there is no condition that can guarantee independence of the initial conditions. The only 
possibility to have independence of the initial conditions is represented by a tauon N 2 - 
dominated scenario [SB], i-e. a scenario where the asymmetry is dominantly produced 
from the next-to-lightest RH neutrinos, and therefore M 2 3> 10 9 GeV, in the tauon 
flavour. The condition M\ <C 10 9 GeV is also important to have projection on the 
orthonormal three lepton flavour basis before the lightest RH neutrino washout |62j . 

5. Density matrix formalism 

As explained in the last Section, leptogenesis is sensitive to the temperature range at 
which the asymmetry is produced. This determines whether the lepton quantum states 
produced in RH neutrino decays either remain coherent or undergo decoherence and 
get projected in flavour space before scattering in inverse processes. Moreover, since 
the lepton states produced in heavy neutrino decays differ in general from the lepton 
flavour eigenstates, lepton flavour oscillations can also in principle arise in a similar 
way as neutrino oscillations happen in vacuum or in a medium. In order to treat the 
problem of flavour oscillations and partial loss of decoherence in a consistent way, one has 
to extend the classical Boltzmann framework to account for these intrinsically quantum 
effects. The formalism of the density matrix is appropriate for this purpose [63] . The 
density matrix for leptons with momentum p is defined as 



with a (a^) denoting the annihilation (creation) operator, i.e. it is the expectation value 
(to be understood in statistical terms) of the generalised number operator. For anti- 
leptons one has analogously a density matrix pi(p). The diagonal elements of the density 
matrix contains nothing else than the occupation numbers of the two flavoured leptons, 
and the off-diagonal elements encode flavour correlations. 




(29) 
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Let us consider an ^-dominated scenario for simplicity. Moreover, let us consider 
a momentum integrated description, introducing the matrix of lepton number densities 
Ne/R 3 = gi J -SSs pi(p), where R is the scale factor. Within the density matrix 
formalism, the asymmetry can be calculated from the following density matrix equation 
in the flavour space t-t l pEU El E3 SHI EQ] 

dN B a ~ L 



dz 



e§D 1 (N Nl -N%] 



-W, {V° m ,N B - L } 



a/3 



+ 1 



. Re(A 7 



Hz 

Im(A T ) 
Hz 






:) 


,N e+! 











B-L 



(30) 



a/3 



where we defined iV a l 



and where V ^ is a matrix projecting the lepton 
quantum states along the flavour '1' of a lepton l\ produced from the decay of a 
RH neutrino N\ . The CP asymmetry matrix is a straightforward generalization of 
Eq. ( )28|) [4*9| |4"8| [60] . The real and imaginary parts of the tau- lepton self-energy are 
respectively given by 



Re(A T ) = &T 



and 



Im(A T ) = 8 x 10 _ Vr r 



(31) 



where f T is the tauon Yukawa coupling. The commutator structure in the third term 
on the RHS of Eq. (1301 accounts for oscillations in flavour space driven by the real part 
of the self energy, and the double commutator accounts for damping of the off-diagonal 
terms driven by the imaginary part of the self energy. 

In order to close the system of equations, we also need an equation for the matrix 
N e+ g, which is given by 



dKf 



Re(A T ) 



where S„ 



l_l . (^UN B - L - S g (Nif - 2 N?6af,) 
Tg/ (Hz) accounts for gauge interactions. As shown in 



(32) 

, this term has 

the effect of damping the flavour oscillations. This can be understood by noticing that 
gauge interactions force Nt = 2N^5 a /s, which in turn renders the oscillatory term 
Eq. (J30J negligible. 

Eq. ( l30i) should be solved in any of the intermediate regimes where lepton states 
are partially coherent. Actually the range of relevance of the density matrix equation 
might be more important than previously believed. Indeed, it was found in [48J that in 
order to recover the unflavoured regime, one should have masses well above 10 13 GeV. 

As already discussed, the contribution from heavier RH neutrinos cannot be 
neglected in general. Therefore, the density matrix equation (130]) should be extended to 
account for such effects; this was done in [60j, where a general equation was presented, 
valid for any RH neutrino decays and any temperature range. Flavour projection effects 
as well as phantom terms are readily taken into account in this framework. 
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6. Limit of quasi-degenerate heavy neutrinos 

If 5 2 <C 1, the CP asymmetries £1,2 get resonantly enhanced as £1,2 oc l/5 2 [3H E2J 133] . 
If, more stringently, 5 2 < 1CT 2 , then r\ B oc l/<5 2 and the degenerate limit is obtained 
[40] . In this limit the lower bounds on Mi and on T re h get relaxed proportionally to S 2 
and at the resonance they completely disappear [39, 67J. The upper bound on mi also 
disappears in this extreme case. In a more realistic case where the degeneracy of the 
RH neutrino masses is comparable to the degeneracy of the light neutrino masses, as 
typically occurring in models with flavour symmetries, the upper bound nil < 0.1 eV 
obtained in the hierarchical case imposing the validity of Boltzmann equations, gets 
relaxed to mi < 0.4 eV [381 EO] in the case M 3 3> M 2 , while it is basically unchanged 
if M 3 = M 2 [30]. The difference is due to the fact that the relaxation is mainly to be 
ascribed to the extra-term in the CP asymmetry mentioned at the end of Section 3. This 
extra-term, subdominant in the hierarchical case, can become dominant (if M 3 3> M 2 ) 
in the quasi-degenerate case (5 2 = (M 2 — M\)jM\ <C 1) and it grows with the absolute 
neutrino mass scale instead of being suppressed as the usual term [38J. On the other 
hand this extra-terms vanishes exactly when M 2 = M 3 , a more reasonable assumption 
for 5 2 < 1. 

In the full three-flavour regime, the contributions from all quasi-degenerate RH 
neutrinos should be taken into account and in this case the final asymmetry can be 
calculated as 

i,a \ j / 

Notice that, for each lepton flavour a, the washout in the degenerate limit is described 
by the sum of the flavoured decay parameters for each RH neutrino species. 

The simplest way to obtain a quasi-degenerate RH neutrino mass spectrum is to 
postulate the existence of a slightly broken lepton number symmetry in the lepton 
sector [681 ESI [70]. Assuming the existence of only two RH neutrinos at first for 
simplicity, it is possible to write down Yukawa couplings and a Majorana mass term 
in a way that conserves lepton number. In the "flavour" basis, which we denote by a 
prime, one RH neutrino can be assigned lepton number +1, and the other one -1, so 
that the seesaw mass matrix from Eq. ([7]) takes the form 



( h' al v 2 \ 
ti£v 2 M 



(34) 

y M ) 
which conserves lepton number. Rotating to the RH neutrino mass basis, one finds that 
/ h' al v 2 ih' al v 2 \ 



My 



hZv 2 M 



(35) 



\ ih'J lV 2 My 

It can be seen that the two RH neutrinos are exactly degenerate in this limit, and one 
can easily show that the neutrino mass matrix m v in Eq. ((8]) vanishes identically. In 
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other words, having non-zero neutrino masses requires a small breaking of the lepton 
number symmetry, which automatically splits the two RH neutrinos into a quasi-Dirac 
fermion pair. 

There are different ways to implement the breaking of the lepton number symmetry, 
thus generating non-zero neutrino masses. For instance, in the above two-RH neutrino 
model, we can write: 



proportional to the breaking parameters e and /j,2, as expected. If the lepton number 
symmetry is broken as in the first term in Eq. (J37l) . it is referred to as 'linear' [73]; if 
the second term is at work, it is referred to as inverse/double seesaw mechanism [54]. 
In the following, we will refer to these models for simplicity as 'inverse seesaw models'. 
No matter how the lepton number symmetry is broken, the bottom line is that these 
models fall in the category of low-scale seesaw models, where the size of the Yukawa 
couplings is not necessarily suppressed if the RH neutrino mass scale is lowered to 
the electroweak scale. This can lead to interesting non-unitarity effects in neutrino 
oscillation experiments [74j [75] , as well as observable lepton flavour violating rates in 
experiments looking for /i — > ej, r — > fi(e)j, \i —> eee or /i -> e conversion in nuclei [44] . 
Note that, within the orthogonal parameterisation Eq. ( fit)]) , the inverse seesaw model 
with two RH neutrinos is obtained in the limit |0| — > oo [76J. 

The simple model in Eq. fl34| can be trivially extended to have a third massive RH 
neutrino of mass still conserving lepton number. It would have zero lepton number and 
be decoupled from leptons. Without further assumptions, the other RH neutrino mass 
scale is independent of M, and therefore it can be much lower or much higher. 

Leptogenesis in the context of these low-scale seesaw models was intensely studied in 
recent years. It was found in [44] that it is possible to have at the same time successful 
leptogenesis and low energy observable effects (beyond standard neutrino oscillation 
phenomenology) such as, for example, charged lepton flavour violation processes. This 
is possible in a model with three RH neutrinos and with the help of very large flavour 
effects (hence the name 'resonant r- leptogenesis'). 

More recently, this possibility was examined in the context of a two-RH neutrino 
model, and it was found that, in the limit M 2 — Mi 3> i\ 2 , leptogenesis does not allow 
large enough Yukawa couplings to have non-trivial consequences at low energies [76J. 
This conclusion was re-examined in [72] relaxing the requirement on the mass splitting, 
and allowing for more extreme quasi-degeneracies. In the limit M 2 — Mi <C r 12 , it was 
argued that the decay parameter K should be replaced by an effective decay parameter 

oc K a {M 2 — Mi) 2 /Y\ which depends explicitly on the small breaking of the lepton 



/ h' al v 2 e a v 2 \ 
M v = h'^v 2 /ii M 

\ t T y m n 2 ) 



(36) 



which implies that the full light neutrino mass matrix is given by [7T| [72] 



m v ~v 2 ( e-U' T + h'j-e T ) - v 2 ( h'j-^j-h'A 
V M M J V M M J 



(37) 
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number symmetry. As a matter of fact, it is expected that the washout of lepton number 
vanishes in the limit of lepton number conservation, and in [72] it was rigorously derived 
from the negative interference between the two RH neutrinos exchanged in the AL = 2 
process ->■ 

A more controversial issue is the behavior of the CP asymmetry parameter in the 
limit M 2 — Mi <C r 12 , which is directly related to the form of the regulator for the RH 
neutrino propagator in the self-energy diagram. The reason is that the location of the 
pole for the RH neutrino determines the maximum enhancement of the CP asymmetry. 
Following [77], Ref. [7J uses 

x (M?-M?)2 + (M i r i -M j r 3 -) 2 (38) 

with the regulator given by the difference MiTi — MjTj, whereas [44J finds the regulator 
MjTj. Within the inverse seesaw model considered (with two RH neutrinos) the decay 
rates I\ and r 2 are predicted to be equal in the lepton number conserving limit. 
Therefore, a regulator M{Ti — M 2 r 2 allows for a much larger enhancement of the CP 
asymmetry than MiTi, implying that leptogenesis is compatible with observable lepton 
flavour violation rates [72]. The precise value of the CP asymmetry in the regime 
M 2 — M\ <C r 12 , which is especially relevant in inverse seesaw models, is currently 
still an open issue, whose resolution presumably lies beyond the classical Boltzmann 
approach (see next Section). 

Note that the Weinberg-Nanopoulos [78] requirement, that at least two couplings 
should violate lepton (or baryon) number to have a generation of asymmetry, is satisfied 
with both regulators in the inverse seesaw model. Indeed, when four out of the five 
lepton-number- violating couplings in Eq. (1361) are turned off, the CP asymmetry vanishes 
with both regulators, as it should [72]. However, the limit of all couplings taken 
simultaneously to zero is not well-behaved for the regulator in Eq. ( 1381) . as noted in 

A different scenario was considered in [55] but still within the inverse seesaw 
framework. There, a third almost decoupled RH neutrino was added, with a mass 
Mi <C M. Lepton number violation was included in the coupling of the lightest RH 
neutrino, such that it can decay producing a lepton asymmetry. Therefore, in this 
case, the quasi-degenerate pair of RH neutrinos is not responsible for the generation 
of asymmetry. Nonetheless, their large couplings to leptons (leading for instance to 
non-unitarity effects in neutrino oscillations) imply that the flavoured CP asymmetry, 
more precisely the lepton number conserving part in Eq. ( l28l) . can be large even 
for TeV-scale RH neutrino masses. However, this scenario of 'non-unitarity driven 
leptogenesis' has intrinsically large lepton flavour violating interactions that lead to 
flavour equilibration [HD]- Using the cross-sections for flavour violating interactions 
obtained in [44], one finds that the asymmetry cannot be generated in the right amount if 
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-^1 ^ 10 8 GeV [53] . However, these cross-sections were found recently to be significantly 
more suppressed than in jH], leading to an interesting relaxation of the bounds down 



It is worth mentioning that leptogenesis was investigated in the context of the 
inverse seesaw model also when lepton number is exactly conserved [81]. In this 
case, the observed baryon asymmetry can be generated if leptogenesis occurs during 
the electroweak phase transition, when the sphaleron rate progressively goes out of 
equilibrium. The lepton flavour asymmetries, generated exclusively thanks to flavour 
effects (again the second term in Eq. (128!) ). are then converted into a baryon asymmetry 
before total washout. 

Another example where leptogenesis occurs in the resonant regime is radiative 
leptogenesis [82l [83]. There, RH neutrinos are assumed exactly degenerate at some 
high scale (for instance the GUT scale), and small RH neutrino mass splittings are 
generated by the running of the seesaw parameters. 

In these low-scale seesaw models with TeV-scale RH neutrinos, it is natural to 
wonder whether some interesting signatures could be observed at collider experiments 
such as the LHC Unfortunately, it seems that within the simplest type-I seesaw, 
the prospects are rather dim [70l [M]. The main problem is that the mixing of 
RH neutrinos with light neutrinos has strong upper limits from rare (lepton-flavour- 
violating) decays [71], which prevents an important production of RH neutrinos at 
the LHC. However, in extended models such as Type II seesaw, Type III seesaw, left- 
right symmetric models, or simply an extra U(1)b-l, the prospects are much more 
encouraging (see [21] in this Issue). 

7. Improved kinetic description 

We have already discussed the density matrix formalism, which goes beyond the 
traditional kinetic treatment with Boltzmann (rate) equations. This Section is devoted 
to other kinetic effects that are important both for an estimation of the theoretical 
uncertainties in the calculation and for a better conceptual understanding of the minimal 
leptogenesis framework. 

7.1. Momentum dependence 

Within the vanilla scenario, the final asymmetry is computed by solving classical 
Boltzmann equations for the RH neutrino and lepton number densities, so-called rate 
equations. These are obtained from the Boltzmann equations for the distribution 
function integrating over momenta with some approximations (see below). One can 
then wonder what is the theoretical error introduced by this integrated description. 

Given a particle species X, the number density is obtained by integrating the 
distribution function over momentum, 



to Mi > 10 6 GeV P 




(39) 
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where gx is the number of degrees of freedom of particle X. For leptogenesis with 
decays and inverse decays, the system of Boltzmann equations (one for the RH neutrino 
and one for lepton number) in the expanding Friedmann-Robertson- Walker Universe is 
given by 

CdUn] (40) 
Cd [ft-e], (41) 



df N 


Pat \n „, 

o\p N 


dt 


dfe-i 


i I u d fe-i 
dpi 


dt 



where the collision integrals on the right-hand side are defined as 

x [f B f c (l-f A )\M(BC^A)\ 2 - (42) 
f A (l-f B )(l-f c )\M{A^BC)\ 2 ] . 

Note that the double-counting problem is solved here in the same way as for the 
integrated Boltzmann equations, by consistently including the resonant part of the 
AL = 2 scatterings. In order to recover the usual Boltzmann equations ( fT6l) -(fl7l). 
one has to introduce three approximations: (i) kinetic equilibrium for the RH neutrinos, 
which can be expressed as /v/Zat 1 — n N/ n< N'i (ii) Maxwell-Boltmann distributions for RH 
neutrinos, leptons and Higgs fields; (iii) neglect Pauli blocking and Bose enhancement 
factors. 

RH neutrinos are only coupled to the thermal bath via their Yukawa couplings. It 
is therefore clear that in the weak washout regime, K\ <C 1, the assumption of kinetic 
equilibrium is not a very good one, and, indeed, it was found that the lepton asymmetry 
computed with the above equations can differ by up to 50% in the weak washout regime 
compared to the usual treatment with integrated equations [851 ESI E7]. However, as 
expected, in the strong washout regime, the above-mentioned approximations are very 
good and the integrated rate equations can be used safely. 



7.2. Non-equilibrium formalism 

Leptogenesis is an intrinsic non-equilibrium problem. One of the Sakharov's conditions 
is indeed that a departure from thermal equilibrium is necessary to produce the 
baryon asymmetry. Therefore, it does not come as a surprise that recently a huge 
effort [HHl [891 [901 EH 1921 [93l SHI [9^1 [95] was made to understand leptogenesis within 
non-equilibrium quantum field theory, also known as the closed-time-path (CTP) or 
Keldysh-Schwinger formalism. This more rigorous, though far more complex, approach 
has the advantage of taking into account quantum effects that are completely missed 
by the usual approach, like memory effects and off-shell effects. Moreover, it allows 
the straightforward inclusion of flavour oscillations and decoherence [48J, it has the 
advantage of being able to consistently account for finite density corrections, and it has 
no double counting problem. 
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Concretely, in the non-equilibrium framework, one needs to find the equations of 
motion for the two-point correlation functions (i.e. the propagators or Green's functions) 
of RH neutrinos and leptons from the general Schwinger-Dyson equation on the CTP, 

Si l (x,y) = S^(x,y) - T, e (x,y) , (43) 
S~\x, y) = S^ix, y) - E N (x, y) , (44) 

which are obtained from the variational principle on the effective action, a functional of 
the full propagators A^, Sg and S, for the Higgs, lepton and RH neutrino, respectively. 
In the above equations, the subscript denotes the free propagators, and and T, N 
are the self-energies for the leptons and RH neutrinos, respectively. Note that the self- 
energies are themselves functions of the propagators. For instance, the one-loop lepton 
self-energy depends on the RH neutrino and Higgs propagators: 

Ef {x, y) = -h ai h^P R S ij (x, y)PiA*{y, x) . (45) 

It usually proves convenient to decompose any two-point function D(x, y) into a spectral, 
D p , and a statistical component, Dp: 

D(x, y) = D F (x, y) - % -sgn(x° - y°)D p (x, y) . (46) 

Convoluting the Schwinger-Dyson equations with the full propagator, we finally arrive 
at a system of two coupled integro-differential equations, the so-called Kadanoff-Baym 
equations: 

i$ x Sip P {x,y) = / d 4 zY, e p 1 (x,z)S(i F B {z,y) 
Jo 

- d 4 zJ: eF \x,z)S i f(z,y), (47) 
Jo 

x S t f(x,y)= / d 4 zJ: e a ;'(x,z)S e f(z,y). (48) 

The corresponding equations for the RH neutrino two-point function are obtained by 
changing $ x — > $ x — M, and S a/3 — > , where a, (3 = e,fi,r are lepton flavours, and 
i,j = 1,2,3 are RH neutrino flavours. It can be noticed that the Kadanoff-Baym 
equations contain an integration over the entire history of the system, a 'memory' 
integral, which encodes all previous interactions with momentum and spin correlations. 
An attempt of studying memory effects in the context of leptogenesis [89] has found 
that large effects could arise in the resonant limit in the weak washout regime. In order 
to recover a Markovian description of the system, characterized by uncorrelated initial 
states at every timestep, one has to perform a gradient expansion (for an alternative 
approach, see [91]), relying on the fact that the microscopic timescale t m i C ~ 1/Mj is 
much smaller than the macroscopic timescales t mac ~ 1/Tj, \/H. This is also known as 
the molecular chaos approximation. 

In leptogenesis, one needs to compute the evolution of the lepton number density. 
The latter is given by the average expectation value of the zeroth component of the 
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lepton number current, given by 

3Lp(z) = -trkfS iafl {z,x)] . (49) 
One then obtains for the lepton number density 

n Laf t(t)=i J (0s jf ^ jf dt " tT [Xt p f(t',t")S eFp (t",t') 

-E lF f(t'X)S ipp (f,t')] , (50) 

after switching to momentum space. From this master equation, one then needs to 
input the equilibrium expression for the lepton and Higgs propagators, as well as the 
non-equilibrium Majorana neutrino propagator. In order to obtain a Boltzmann-like 
equation, two further simplifications are required. First, the quasi-particle ansatz, also 
known as the on-shell approximation, which states: 

D p {X,p) = 2vrsgn(p°)5(p 2 - m 2 ) , (51) 

where X = (x + y)/2 is the central coordinate. Then, one can express the statistical 
propagator for the occupation number using the so-called Kadanoff-Baym ansatz, 

r 



D F (X,p) 



f(x, P ) + 2 



DJX,p), (52) 



which is chosen such that, in equilibrium, the correct fluctuation dissipation relation 
between D e p{p) and D^(p) is automatically obtained (see [90] for more details). With 
these approximations, one arrives at a Boltzmann equation which includes finite density 
effects. This equation would be similar to the kinetic equation ( I4ip . with however the 
important new property: 

\M{N ^l§)\ 2 = |M(^> ^ iV)| 2 = \M (N £$)| 2 (1 +£(p,T)), (53) 

where subscript '0' denotes tree level, and [92] 

e(p,T) =ex(l + J ^[f; q m - fP(E 2 )fj (54) 

where e is the CP asymmetry defined in Eq. (fl4l) . and Eip = | (M 2 + p 2 )^ ±pcos9 . 
Note that the above property in Eq. (|53|) explicitly avoids the double counting problem, 
which plagues the momentum-dependent description of the last Subsection (as well 
as the usual integrated one). As it can be seen, the CP asymmetry parameter now 
includes finite density effects via a dependence on the Higgs and lepton distribution 
functions. This result agrees with that one found using thermal field theory in the real 
time formalism, when the right convention is used [92J, as well as in the imaginary 
time formalism as recently obtained [96J. It, however, disagrees with an earlier work 
based on thermal field theory in the real time formalism, where a term quadratic in the 
distribution functions was found [29] . 

One should then wonder whether the non-equilibrium formalism described above 
leads to substantial differences in the leptogenesis predictions for the baryon asymmetry. 
It turns out that the modifications are very important (up to an order of magnitude) 
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in the regime where Mj/T <C 1, but they tend to vanish in the non-relativistic 
regime Mj/T ^> 1. In other words, major changes are expected mainly in the weak 
washout regime for a vanishing initial RH neutrino abundance, because in this case 
the asymmetry is produced in both regimes: the initial asymmetry is produced with the 
wrong sign at early times and is compensated by the right-sign asymmetry at later times. 
This was confirmed in [93], where it was found that the sign of the final asymmetry 
could even get changed by finite density effects. In all other cases, corrections are 
quantitatively small, in particular in the strong washout regime, Ki 3> 1, where the 
asymmetry is produced exclusively in the non-relativistic regime [92j [93]. 

7.2.1. Resonant limit It is interesting to study whether the above quantum kinetic 
formalism can provide some insight onto the CP asymmetry produced in the extreme 
quasi-degenerate limit, \M\ — M 2 | <C M lj2 . This formalism includes off-shell effects, 
as well as the proper inclusion of coherent transitions Ni — > Nj, both of which 
should be important when \M\ — M 2 | ~ I\ 2 . This problem attracted some attention 
recently [941195]. 

In [91] it was found that the enhancement of the CP asymmetry derived within the 
Kadanoff-Baym (KB) formalism, 

pKB = M x M 2 {Ml - Ml) 

~ (Ml-Mfy + (M 1 T 1 + M 2 T 2 ) 2} { ' 

differs from the Boltzmann result, 

.be = M 1 M 2 (M* - Ml) 

~ (Mf - Mlf + {M 1 Y 1 - M 2 r 2 )2 ' {00) 

specifically because of a contribution from coherent RH neutrino oscillations. This would 
prevent any additional enhancement when Ti ~ T 2 , as predicted by the inverse seesaw 
model. In such a case, the new region in the parameter space of successful leptogenesis 
explored in [72] is not available. 



7.2.2. Adding flavour We saw in Section 5 how the first kind of 'quantum' effects 
were included, namely flavour oscillations and decoherence, by switching from classical 
Boltzmann equations for lepton number densities to evolution equations for the density 
matrix. For the lightest RH neutrino and two relevant lepton flavours, we presented 
Eq. ( 1301 . The CTP formalism can also account for a flavour matrix structure, through 
the lepton propagators and self-energies. Reassuringly, the structure of Eq. fl3TJ|) was 
also found within this formalism [4"8] . 

8. Other corrections 

8.1. Thermal effects 

Leptogenesis occurs in the very hot thermal bath of the early Universe. Leptons 
and Higgs fields have very fast interactions with the thermal bath due to their gauge 
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couplings. This gives them an effective mass which is proportional to temperature |65j. 
This effective mass allows processes which were otherwise kinematically forbidden to 
occur, such as CP-violating Higgs decays to RH neutrinos, when T 3> M. Thermal 
effects have therefore a direct impact on leptogenesis, and the first study to try and 
quantify them (29] employed a real time formalism and hard thermal loop resummation. 
It was found that the CP violating parameter has a strong temperature dependence, and 
that in the weak washout regime, there could be important differences with the zero- 
temperature treatment. On the other hand, in the strong washout regime, if>l, the 
usual results from a vacuum calculation are recovered with a good accuracy. Recently, 
a new study came out using the imaginary time formalism [96] , and the CP asymmetry 
parameter was found to differ from (29], and to agree with more recent attempts with 
non-equilibrium quantum field theory (see previous Section). 

8. 2. Spectator processes 

Chemical equilibrium holds among Standard Model particles in the early Universe 
above the electroweak phase transition thanks to gauge interactions. On the other 
hand, Yukawa interactions only force some new conditions (between left- and right- 
handed fermions) when the temperature is high enough, depending on the size of the 
Yukawa coupling. If one imposes additionally hypercharge neutrality and the effect of 
electroweak and strong sphaleron equilibrium, one arrives at a set of relations among 
the chemical potentials (or among the asymmetries) of leptons, Higgs and baryons [97J. 
These processes, although not directly involved in the leptogenesis process, hence the 
name 'spectator processes' [98], have an indirect effect on the final asymmetry through 
a modified washout. The first effect is the inclusion of the Higgs asymmetry as a 
new contribution to the washout. The second effect is that the asymmetry originally 
produced in leptons of a particular flavour gets redistributed into the other flavours 
following precise relations [26] . This means that the set of Boltzmann equations to solve 
for the different lepton flavours are now coupled to each other via a flavour coupling 
matrix. 

Within the iVi-dominated scenario, the overall effect of spectator processes is 
usually subdominant compared to flavour effects, but it can change the final result by as 
much as 40% depending on the temperature at which the asymmetry is produced [971157]. 
Supersymmetry has new degrees of freedom and new constraints can be derived. A full 
study was performed recently in |99j and the overall effect was found to be again of 
order one. Within the ^-dominated scenario, however, potentially much bigger effects 
are possible [59J. 

8. 3. Scattering processes 

It can be shown that leptogenesis is well-described in the strong washout regime, K 1, 
by just decays and inverse decays. The reason is that in this regime the asymmetry 
is produced at relatively late times, at T <C M, with no dependence on the dynamics 
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happening at T > M when scattering processes are important. However, they should be 
included in the weak washout regime for initial vanishing RH neutrino abundance. For 
instance, AL = 1 Higgs-mediated scatterings involving top quarks, such as £N -h- Q^i, 
contribute to the washout and the CP asymmetry [19], and their inclusion is crucial to 
have a correct estimation of the final asymmetry in the weak washout regime |100j . 
Scatterings involving gauge bosons, such as £N <H- have also been included, 

especially their contribution to the CP asymmetric source term jlOlj . It was found 
that the factorization of the CP asymmetry from decays and scatterings involving top 
quarks does not happen with scatterings involving gauge bosons. In particular, there is 
a new source of lepton-number-conserving CP asymmetry. 

9. Testing new physics with leptogenesis 

The seesaw mechanism with three RH neutrinos extends the Standard Model by 
introducing eighteen new parameters. On the other hand, low-energy neutrino 
experiments can only potentially test the nine parameters in the low-energy neutrino 
mass matrix m u . Nine high-energy parameters, those characterising the properties of 
the three RH neutrinos, the three masses and the six parameters encoded in the seesaw 
orthogonal matrix of Eq. (flOl) . basically fixing, together with the light neutrino masses, 
the three lifetimes and the three total CP asymmetries, are not tested by low-energy 
neutrino experiments. Quite interestingly, the requirement of successful leptogenesis, 



provides an additional constraint on a combination of both low-energy neutrino 
parameters and high-energy neutrino parameters. However, just one additional 
constraint would not seem sufficient to over-constrain the parameter space leading 
to testable predictions. In spite of this observation, as we have seen, in the vanilla 
leptogenesis scenario one can derive an upper bound on the neutrino masses. The reason 
is that, within this scenario, the dependence of t\b on the six parameters related to the 
properties of the two heavier RH neutrinos cancels out. In this way the asymmetry 
depends on a reduced subset of high-energy parameters (just three instead of nine). At 
the same time, the final asymmetry gets strongly suppressed when the absolute neutrino 
mass scale it is larger than the atmospheric neutrino mass scale. For all these reasons, 
by maximising the final asymmetry over the high energy parameters and by imposing 
successful leptogenesis, an upper bound on the neutrino masses is found. 

When flavour effects are considered, the vanilla leptogenesis scenario holds only 
under very special conditions, as we have seen. In general, the final asymmetry depends 
also on the parameters in the leptonic mixing matrix. Therefore, accounting for flavour 
effects, one could naively hope to derive definite predictions on the leptonic mixing 
matrix as well, in addition to the upper bound on the absolute neutrino mass scale. 
However, the situation is quite different when flavour effects are taken into account. 
This is because the final asymmetry depends, in general, also on the six parameters 




(57) 
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related describing the two heavier RH neutrino properties and that were cancelling out 
in the calculation of the final asymmetry in the vanilla scenario, and this goes at the 
expense of predictability. 

For this reason, in a general scenario with three RH neutrinos and flavour effects 
included, it is not possible to derive any prediction on low-energy neutrino parameters. 
As we discussed, even whether the upper bound rrii < 0.1 eV on neutrino masses still 
holds or not is an open issue and a precise value seems to depend on a precise account 
of many different subtle effects and in particular it necessarily requires a density matrix 
formalism. 

In order to gain predictive power, two possibilities have been explored in the past 
years. 

A first possibility is to consider non-minimal scenarios giving rise to additional 
phenomenological constraints. For example, as discussed in Section 6, the inverse seesaw 
model, which technically can be still regarded as part of the minimal type I seesaw 
model, allows for non-trivial phenomenologies at low energy beyond standard neutrino 
oscillations, such as non-unitarity effects and observable lepton flavour violation. An 
experimental observation of any of these signatures would be of great value to reduce the 
freedom in the choice of seesaw parameters. In recent years, during the Large Hadron 
Collider era, it has been also intensively explored the possibility that, within a non- 
minimal version of the seesaw mechanism, one can have successful low scale leptogenesis 
together with collider signatures. It has also been noticed that in the supersymmetric 
version of the seesaw, the branching ratios of lepton-flavour- violating processes or electric 
dipole moments are typically enhanced, and hence the existing experimental bounds 
further constrain the seesaw parameter space [1021 1103] . 

A second possibility is to search for a reasonable scenario where the final asymmetry 
depends only on a reduced set of independent parameters over-constrained by the 
successful leptogenesis condition, as with the vanilla scenario. From this point of view, 
the account of flavour effects has opened very interesting new opportunities or even re- 
opened old attempts that fail within a strict unflavoured scenario. Let us briefly discuss 
some of the main ideas that have been proposed within this second possibility, the first 
one being covered elsewhere in this Issue [2~T] . 

9.1. Two-RH neutrino model 

A phenomenological possibility that attracted great attention is the two-RH neutrino 
model [104J, where the third RH neutrino is either absent or effectively decoupled in 
the seesaw formula. This necessarily happens when M 3 3> 10 14 GeV, implying that the 
lightest left-handed neutrino mass mj has to vanish. It can be shown that the number 
of parameters decreases from 18 to 11 in this case. In particular the orthogonal matrix 
is parameterised in terms of just one complex angle. 

In leptogenesis the two-RH neutrino model has been traditionally considered as a 
sort of benchmark case for the iVi-dominated scenario, where the final asymmetry is 
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dominated by the contribution from the lightest RH neutrinos [1HJ EH 1105] . However, 
recently, as we anticipated already, it has been shown that there are some regions in the 
one complex angle parameter space that are ^-dominated |106] and that correspond 
to so-called light sequential dominated models [53] . 

It should be said that even though the number of parameters is highly reduced, in 
a general two-RH neutrino model it is still not possible to make predictions on the low- 
energy neutrino parameters. To this extent, one should further reduce the parameter 
space, for example assuming texture zeros in the neutrino Dirac mass matrix. 

9.2. SO (10) -inspired models 

In order to gain predictive power, one can impose conditions within some model of 
new physics embedding the seesaw mechanism. An interesting example is represented 
by the 'S'O(10)-inspired leptogenesis scenario' |107[ I108j 1109] . where SO(10)-inspired 
conditions are imposed on the neutrino Dirac mass matrix mo- In the basis where 
the charged leptons mass matrix and the Majorana mass matrix are diagonal, in the 
bi-unitary parametrisation, one has mo = Vl D mD Un, where D mD = diag(Ai, A2, A3) 
is the diagonalised neutrino Dirac mass matrix and the mixing angles in Vl are of 
the order of the mixing angles in the CKM matrix Vckm- The Ur and three Mj 
can then be calculated from V^, U and m^, since the seesaw formula Eq. (jHJ) directly 
leads to the Takagi factorisation of M~ l = D^ V^U D m U T Vl D^}, or explicitly 
M-^ = U R D- A }Ul 

In this way the RH neutrino masses and the matrix Ur are expressed in terms of 
the low-energy neutrino parameters, of the eigenvalues A« and of the parameters in Vl- 
Typically one obtains a very hierarchical spectrum Mi ~ 10 5 GeV and M 2 ~ 10 11 GeV, 
the asymmetry produced from the lightest RH neutrino decays is by far unable to explain 
the observed asymmetry [108] . 

However, when the N 2 produced asymmetry is taken into account, successful (N 2 - 
dominated) leptogenesis can be attained [110] . In this case, imposing the leptogenesis 
bound and considering that the final asymmetry does not depend on Ai and on A 3 , one 
obtains constraints on all low-energy neutrino parameters, which have some dependence 
on the parameter A2 typically parameterised in terms of a 2 = ^/^c, where m c is 
the charm quark mass. Some examples of the constraints on the low-energy neutrino 
parameters are shown in Fig. 3. They have been obtained scanning over the 2a ranges of 
the allowed values of the low-energy parameters and over the parameters in Vl assumed 
to be / < Vl < Vckm and for three values of a 2 = 5,4, 1. It is particularly interesting 
that when the independence of the initial conditions is imposed, negative values of Jcp 
seem to be favoured [112] establishing a connection between the sign of Jcp and of the 
matter-antimatter asymmetry. 

A supersymmetric version of this scenario including the renormalization group 
evolution of all the relevant couplings was also studied in [113] . and including a type 
II contribution to the seesaw mechanism from a triplet Higgs in left-right symmetric 
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models in |114j . 

Together with strongly hierarchical RH neutrino mass patterns, there exist also 
'level crossing' regions where at least two RH neutrino masses are arbitrarily close to each 
other, in a way that the lightest RH neutrino mass is uplifted and CP asymmetries are 
resonantly enhanced to some lever. Also in this case, imposing the successful leptogenesis 
condition, one then obtains conditions, though quite fine tuned ones, on the low energy 
neutrino parameters |109j . 

9. 3. Discrete flavour symmetries 

An account of heavy neutrino flavour effects is also important when leptogenesis is 
embedded within theories that try to explain tribimaximal mixing for the leptonic 
mixing matrix via flavour symmetries. It has been shown in particular that, if the 
symmetry is unbroken, then the CP asymmetries of the RH neutrinos would exactly 
vanish. On the other hand, when the symmetry is broken, for the naturally expected 
values of the symmetry breaking parameters, then the observed matter-antimatter 
asymmetry can be successfully reproduced |115[I116[I117[[TT8] . It is interesting that in a 
minimal picture based on an A4 symmetry, one has a RH neutrino mass spectrum with 
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10 15 GeV > M 3 > M 2 > M l > 10 12 GeV. One has therefore that all the asymmetry 
is produced in the unflavoured regime and that the mass spectrum is only mildly 
hierarchical (it has actually the same kind of hierarchy as light neutrinos). At the 
same time, the small symmetry breaking imposes a quasi-orthogonality of the three 
lepton quantum states produced in the RH neutrino decays. Under these conditions 
the washout of the asymmetry produced by one RH neutrino species from the inverse 
decays of a lighter RH neutrino species is essentially negligible. The final asymmetry 
then receives a non-negligible contribution from the decays of all three RH neutrinos 
species. 

9.4- Super symmetric models 

Within a supersymmetric vanilla framework the final asymmetry is only slightly modified 
compared to the non-supersymmetric calculation [119J. However, supersymmetry 
introduces a conceptually important issue: the stringent lower bound on the reheat 
temperature, T reh > 10 9 GeV, is typically marginally compatible with an upper bound 
from the avoidance of the gravitino problem T re h < 10 6-10 GeV, with the exact value 
depending on the parameters of the model [1201 H21j 1122] . It is quite remarkable that 
the solution of such an issue inspired an intense research activity on supersymmetric 
models able to reconcile minimal leptogenesis and the gravitino problem. Of course, on 
the leptogenesis side, some of the discussed extensions beyond the vanilla scenario that 
relax the RH neutrino mass bounds also relax the T re h lower bound. However, notice 
that in the N 2 dominated scenario, while the lower bound on Mi simply disappears, 
there is still a lower bound on T re h that is even more stringent, T re h > 6 x 10 9 GeV |35j. 

As we mentioned already, with flavour effects one has the possibility to relax the 
lower bound on T re h if a mild hierarchy in the RH neutrino masses is allowed together 
with a mild cancellation in the seesaw formula [30] ■ However for most models, such as 
sequential dominated models [53J, this solution does not work. A major modification 
introduced by supersymmetry is that the critical value of the mass of the decaying RH 
neutrinos setting the transition from an unflavoured regime to a two-flavour regime 
and from a two-flavour regime to a three flavour regime is enhanced by a factor 
tan 2 j3 [SI 1123] . This has a practical relevance in the calculation of the asymmetry within 
supersymmetric models and it is quite interesting that leptogenesis becomes sensitive 
to such a relevant supersymmetric parameter. Recently, a refined analysis, mainly 
discussing how the asymmetry is distributed among all particle species, has shown 
different subtle effects in the calculation of the final asymmetry within supersymmetric 
models finding corrections below 0{1) [99] . 

10. Future prospects for testing leptogenesis 

In 2011 two important experimental results have been announced that, if confirmed, 
can be interpreted as positive for future tests of leptogenesis. 
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The first result is the discovery of a non- vanishing 13 ~ 9° (cf. Section 2) confirmed 
now by various experiments (both long-baseline and reactor). An important consequence 
of the measurement of such a 'large' #13 is the encouraging prospects for the discovery of 
the neutrino mass ordering (either normal or inverted) in current or near-future neutrino 
oscillation experiments such as T2K and NOi^A. For instance, if the ordering is found 
to be inverted in these experiments, we should expect a signal in 0v(3{3 experiments in 
the next decade. If it is not found, the Majorana nature of neutrinos, and therefore the 
seesaw mechanism, would be ruled out. 

Moreover, such a large value of #13 opens the possibility of a measurement of 
the neutrino oscillation CP- violating invariant Jcp oc sin #13 sin 5 during next years. 
This would have some direct model-independent consequences. If a non-vanishing and 
close-to-maximal value is found (|sin5| ~ 1), even the small contribution to the final 
asymmetry uniquely stemming from a Dirac phase could be sufficient to reproduce the 
observed final asymmetry. More generally, the presence of CP violation at low energies 
would certainly support the presence of CP violation at high energies as well, since given 
a generic theoretical model predicting the neutrino Dirac mass matrix mr>, these are in 
general both present. 

A more practical relevance of such a measurement of Jcp, even if in the end it 
indicates a vanishing value within the experimental error, is that it will provide an 
additional constraint on specific models embedding the seesaw, satisfying successful 
leptogenesis and able to make predictions on the low energy neutrino parameters. In 
this way, the expected improvements in low-energy neutrino experiments, made easier 
by a non- vanishing # 13 , will test the models more and more stringently. 

The second important experimental result of 2011 is the hint for the existence of the 
Higgs boson reported by the ATLAS and CMS collaborations. There are at least two 
reasons why this hint can be also interpreted as positive for leptogenesis. First, because 
the whole leptogenesis mechanism relies on the Yukawa coupling between Higgs, lepton 
and RH neutrino. Second, because the measurement of the Higgs boson mass could in 
future open opportunities for additional phenomenological information to be imposed 
on leptogenesis scenarios, for example relying on the requirement of Standard Model 
electroweak vacuum stability [124J. In this respect, it is interesting to note that the 
current value of 125 GeV is compatible with reheating temperatures as high as 10 15 GeV 
as needed by thermal leptogenesis, and it is also compatible with the requirement that 
the Yukawa couplings of the RH neutrinos do not destabilise the Higgs potential when 
RH neutrino masses Mj, assumed to be quasi-degenerate, are smaller than 10 14 GeV. 

What are other possible future experimental developments that could further 
support the idea of leptogenesis? An improved information from absolute neutrino 
mass scale experiments, both on the sum of the neutrino masses from cosmology and 
on m ee from 0uf3f3 experiments could be crucial. For instance, if cosmology provides 
a measurement of the neutrino mass 0.01 eV < m\ < 0.2 eV in the next few years, as 
it is reasonable to expect, then a positive signal in 0vf3f3 experiments must be found, 
otherwise Majorana neutrinos and the seesaw mechanism will be disfavored. On the 
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other hand, in case of a positive signal in experiments, we will be able to say 

that minimal leptogenesis with hierarchical heavy neutrino masses works in a optimal 
neutrino mass window 10 _3 eV < mi < 0.1 eV [2H1 [3QJ EE], where independence of the 
initial conditions is more easily obtained and where we know that successful leptogenesis 
can be safely obtained from existing calculations using Boltzmann equations. Moreover, 
a determination of the allowed region in the plane m ee -^ i could provide an 
additional test of specific leptogenesis scenarios, such as, for example, SO(10)-inspired 
scenarios. 

In conclusion, we are living in an exciting time where new experimental information 
is coming, which is providing and will continue to provide crucial tests for new physics 
models, in particular leptogenesis, in the near future. 
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